clear all; close all; clc;

% % Внутренние параметры
R = 120; % Ом
C = 33e-12; % Ф
L = 100e-6; % Гн
% 
omega0 = 1/sqrt(L*C);
% Резонансная частота
f0 = omega0 / 2 / pi;
fprintf('f0 = %f MHz\n', f0/1e6);
% 
% % Ось времени
T = 1 / (10*f0);
t = 0:T:(15 * 1/f0);
lt = length(t);
% 
% % Выделение памяти
 Uc = nan(1, lt);
 Ur = nan(1, lt);
 dUc = nan(1, lt);
 i = nan(1, lt);
 di = nan(1, lt);
 U = nan(1, length(t));
 dU = nan(1,length(t));
 
% % Внешние параметры, внешние фазовые переменные
f = f0/2:((2*f0-f0/2)/100):2*f0;
A = 1;
K = nan(1, length(f));
for j= 1:length(f)
    E = A*sin(2*pi*f(j)*t);
    
    for k = 2:length(t)
         i(k)=i(k-1) + di(k-1)*T;
        Uc(k) = U(k-1) + dU(k-1)*T;
%         dU(k) = (E(k) - U(k)) / RC;
        Ur(k)= E(k)- Uc(k) - L*i(k);
%         i(k)=i(k-1) + di(k-1)*T;
    end
    
    Us = Ur(fix(end/2):end);
    K(j) = (max(Us) - min(Us))/2;
end

figure;
plot(f,K)
xlabel('f,Hz');
ylabel('K, dB')
grid on
% K = nan(1, length(f));
% for j= 1:length(f)
% E = A*cos(2*pi*f(j)*t);
% 
% % Начальное состояние
% Uc(1) = 0;
% 
% % Начальное напряжение конденсатора
% i(1) = 0;
% 
% % В начальный момент ток в цепи отсутствует
% di(1) = 0;
% Ur(1) = di(1) * R;
% 
% % Напряжение на катушке
% for k = 2:lt
% 
% % % Напряжение на диоде
% % if (E(k-1) - Ul(k-1) - Uc(k-1)) < Uo
% % Ud(k-1) = Rd*i(k-1);
% % else Ud(k-1) = rd*i(k-1) + Uo;
% % end
%  Ur(k) = E(k) - Uc(k-1) - di(k-1)*T*L ;
% 
% % % Напряжение на катушке
% %  di(k) = Ul(k) / L;
% 
% % вызывает приращение тока,
% i(k) = i(k-1) + di(k-1)*T;
% dUc(k) = i(k) / C;  % который определяет заряд конденсатора
% Uc(k) = Uc(k-1) + dUc(k)*T;
% K(k)= Ur(k)/E(k);
% end
% 
% end
% figure;
% plot(f,K(k))
% grid on
% i = sqrt (-1) ; 
% k= (R*sqrt(C/L)*j*2*pi*f)/(1+((j*2*pi*f)^3)*L*C+1+j*2*pi*f*R*C)
% for l=1:length(f)
%     b(l)=abs(k(l));
% end
% plot (f,b)